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Abstract 

We investigate interface motion in disordered ferromagnets by means of Monte Carlo 
simulations. For small temperatures and driving fields a so-called creep regime is 
found and the interface velocity obeys an Arrhenius law. We analyze the corre- 
sponding energy barrier as well as the field and temperature dependence of the 
prefactor. 
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1 Introduction 



Using the random-field Ising model (RFIM) we study the dynamics of driven 
interfaces in disordered ferromagnetic systems by means of Monte Carlo sim- 
ulations. The interface separates regions of opposite magnetization and it is 
driven by a homogeneous field H . Without thermal fluctuations (T = 0) its 
dynamics is affected by the disorder in the way that the interface moves only 
if a sufficiently large driving field is applied. With decreasing driving field the 
interface velocity vanishes at a critical threshold H c and it is pinned below 
H c , due to the random-field. The field dependence of the interface velocity is 
given by 



v{h) ~ hP (1) 
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with h = (H — H c )/H c > 0, i.e. the pinning/depinning transition can be 
considered as a continuous phase transition. 

For finite temperatures T > no pinning of the interface occurs since the 
energy needed to overcome local pinning centers can be provided by thermal 
fluctuations. In particular for driving fields well below H c the interface velocity 
remains finite and is expected to exhibit a so-called creep motion in which 
the interface velocity obeys an Arrhenius law. This Arrhenius law has been 
found for the Edwards- Wilkinson equation with quenched disorder [1,2] and 
we compare these results with those exhibited by a driven interface in the 
3d RFIM. The Hamiltonian of this model is given by 



The sum in the first term is taken over all pairs of neighbored spins. It describes 
the exchange interaction between neighboring spins (Si = ±1) and a parallel 
alignment of nearest neighbors is energetically favored. Additionally, the spins 
are coupled to a homogeneous driving field H and to a quenched local random- 
field hi with (hihj) oc 8^ and {hi) = 0. The probability density p that the 
local random-field takes some value hi is given by p(hi) = (2A) _1 for \hi\ < A 
and zero otherwise. Using simple cubic lattices with antiperiodic boundary 
conditions we started each simulation with an initially flat interface. In non- 
disordered systems the interface moves for any finite driving field. This limiting 
behavior can be recovered if the interface moves along the diagonal direction 
of simple cubic lattices (see for details [3]). In our Monte Carlo simulations 
the dynamics of the interface motion is given by a Glauber dynamics with 
transition probabilities according to a heat-bath algorithm (see e.g. [4]). The 
basic quantity in our investigations is the interface velocity. Since a moving 
interface corresponds to a magnetization M which increases with time t (given 
in Monte Carlo steps per spin) the interface velocity is obtained according to 
v ~ (dM/dt) where (...) denotes an appropriate disorder average. 



2 Creep motion 

The depinning transition occurs at zero temperature only. Including thermal 
fluctuations no pinning of the interface takes place. At the critical threshold 
itself the velocity grows with the temperature according to v (H = H c ) ~ 
T 1 / 5 [3,5]. Well below H c the interface velocity is expected to obey an Arrhe- 
nius law 



H = — — ^2 SiSj - H -Yl^iSi ■ 



(id) 



(2) 



v ~ C(H,T) exp 



E{H) 



(3) 
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Fig. 1. Interface velocities as a function of the inverse temperature for different values 
of the driving field H = 0.3, 0.35, 0.4, 0.6. The velocities are rescaled according 
to Eqs. (3,4). Varying x we obtain nearly straight lines for x = 0.79 ± 0.09. 

which is characterized by its prefactor C (H, T) and an effective energy bar- 
rier E(H). Recent investigations of the creep regime in the QEW yield that 
the prefactor of the Arrhenius law is given by [1] 

C{H,T) = c{H)T- x (4) 

with c(H) being some function of the driving field and the temperature de- 
pendence is characterized by an exponent x. We apply this ansatz to our 
numerically obtained interface velocities in the 3d RFIM. Figure 1 shows the 
data rescaled according to Eqs. (3,4). As can be seen from the data, we ob- 
tain nearly straight lines for x = 0.79 ± 0.09 suggesting that the temperature 
dependence of the prefactor as well as the Arrhenius law itself are a proper 
description of the velocity in the creep regime. The slope of the lines in Fig. 1 
corresponds to the effective energy barrier E(H). Plotting \nvt x vs. E(H)/T 
the different curves collapse onto a single curve (see Fig. 2) without making 
any assumption about c(H). This suggests that c(H) const. 



3 Summary 

We study the creep regime exhibited by a driven interface in the random-field 
Ising model. The creep regime occurs for driving fields well below the critical 
threshold and small but finite temperatures. We find that the interface veloc- 
ity obeys an Arrhenius law. Inspired by a renormalization group approach [1] 
we investigate the details of this Arrhenius law. Studying the corresponding 
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Fig. 2. Rescaled interface velocities. The data coincide with those shown in Fig. 1. 
They are rescaled with the numerically determined energy barrier E(H). The coa- 
lescence of the data occurs without any consideration of c(H). This indicates that 
c(H) k, const. 

prefactor [Eq. (4)] we find that its details do not coincide with [1]. In particu- 
lar, these authors observed a significant field dependence of the corresponding 
prefactor c(H) ~ with ji > and obtained a negative exponent x. Es- 
pecially the opposite sign of x is remarkable and further investigations are 
needed to clarify this point. 
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